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Abstract 

We show that the plasma frequency in wire photonic crystals depends upon the Bloch vector. 
An accurate formula is given. 
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I. INTRODUCTION 
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Photonic crystals are artificial materials structured periodically, with the aim at con- 
trolling light propagation in them Metallic photonic crystals Q] have recently at- 
tracted the attention of several authors for their ability of exhibiting a very low plasma 
requencv and a homogeneous behavior corresponding to a negative relative permittivity 

a negative permeability 
lol . l2o| . The mixture of both media lead to left-handed 
Soukoulis et al. have inspected various formulae for 



materials j2l|. In a recent paper 
deriving the homogenized permittivity and the cut frequency for a system of parallel metal- 
lic wires. Various approaches have been proposed, all of them leading to different formulas, 
none of which being really fully satisfactory. In this letter we provide an accurate formula for 
the cut-frequency as a function of the angle of incidence (or, more generally, the horizontal 
component of the Bloch vector), for a wire mesh photonic crystal. The wire mesh photonic 
crystal is modelized as a stack of diffraction gratings made of very thin rods (radius a) with 
period d (see fig. 1). It is illuminated by an incident, s-polarized, plane wave E'^ under the 
incidence 6. The wavelength in vacuum is A and the wavenumber is = 27t/X. We denote 
by the total electric field and by E'^ = E^ — E"^ the diffracted field. The basic assumption 
here is that koa <^ 1 and k^d <^ 1. 



II. ASYMPTOTIC ANALYSIS OF THE SCATTERED FIELD 

The field diffracted at r = [x, y) by the wire situated at abscissa n x d reads as: E^{r) = 
bnHQ^\kQ\r — nde^l). Here the wires are periodically settled and thanks to Bloch theorem we 
have: 6„ = g^'^'^orfsme^^^ where bo is a coefficient to be determined. Adding the contribution 
of all wires, we get, for the total diffracted field: 

E'^{r) = 6o5^e"^°'^^^'^^i7^^^(A;o|r - nde^\) (1) 
1—1 " 

From scattering theory j2^, we know that the field diffracted by one rod is obtained 
through the scattering matrix So(fcoa) = — Jo{kQa) / HQ^\kQa) by: bn = So[E''{ndex) + 
Y.m^nbmHl{ko\m - n\d)], that is: 

^^(0) 



(2) 



1^0 
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(Jo and Hq^^ are respectively the Bessel function of O*'' order and the Hankel function of 
first type and O*'^ order). 

The total diffracted field is given by: 

(r) = 60 E e*"^°'^^'°^f/i^^ {ko\r~ nrfx|) (3) 

n 

Our aim is now to study the limit of this expression when the wavelength is large with 
respect to the radius of the rods and the period of the grating. The fundamental quantities 
are then k^a and k^d. For later purpose, we define the following combined quantities: 

where /3o = ko cos 9. 

The point is to evaluate the asymptotic behavior of 60 and that of the series defining 
E^(r). For that purpose, we note the following representation formula: 

n " n Pn 



where /3„ = y /cq — a^, a™ = ko sin 6 + n-^. 

The behavior of Bq when k^d tends to cannot be obtained by taking the limit of each 
term of the series in (0), because the terms are singular there. Instead, we perform an 
asymptotic analysis, starting with equality (0)- Our point is to let r = {x,y) tend to 0, and 
we first set y = 0: 

J2 e-"o'^i7(i) (fco \x - nd\) = -Hi'\ko\x\) + W ^e'"-' (6) 

We cannot let x tend to in the right hand side, because HQ^\ko\x\) and the series are 
singular there. However, both singularities compensate. Indeed, the singularity of the series 
can be obtained explicitely by analysing the convergence ratio of the series: as n tends to 
infinity we have /3„ ~ 2z7r|n|/(i, so the series is logarithmic: 

2 1 2 2 / 1 d \ 1 

dnPn df3o d^^ypri 2m\n\ ) ^o^'^H 

and the last series is easily seen to be equal to : ^e*°°^ln[2sin(7rx/(i)]. By using the 
expansion of HQ'\ko\x\) near a; = 0, we get (7 is the Euler constant): 

-Hl,'\ko\x\) + e-o-|ln[2sin(7rx/rf)] ~ -1 - f 7 + f ln(f ) (7) 



which shows that : 

E e--^i/« {ko \n\ d) = -1 - ^7 + - In(^) + ;^ + § E f 7^ + 75^ " (8) 
Finally, when k^d is small, the last series is equivalent to 

and consequently 

E e-^^'H^^ (ko \n\ d) = -1 - -7 + - In A + ^ + O{{kodf) (9) 

We are now in a position to give an asymptotic expansion for 60 • When k^a tends to 
zero, we have: 

s-,\koa) = -^TI^ = -1 - + ln(W2)) + 7^(M' + O{{koa)') . (10) 

From this expression and where we remove the terms that tend to with k^d, we obtain 
from 

5g ^ _^^^^*(0) (11) 

We are now able to give an asymptotic expansion of E'^ by splitting the sum Q into two 
terms corresponding to the evanescent modes and the unique propagative mode : 

E'^ r-. E'^ + E'^ (12) 

prop ' evan \ J 



where 



E^p = e^("°^+^°'^l)i?^(0) , i?Ln = E ^e^(-^+^"'^l)i?^(0) . (13) 



The total field is obtained by adding the incident field: 

For y > : E\r 
For y <0 : E\r 

The reflection and transmission coefficients are readily obtained: 



„i(anx-0r,v) , -1 i(anx+0()v) 1 X J_gi(Q:na;+/3n|j/|) _g''YO) 

X „i(anx-0nv) 1 X , J_p«{ona;+/3„|y|)l E^iffl'l 

X+i L 1-ix ^"7^0 J -C' lu; . 



t = l + r = ^- (15) 



1 - ZX ' ' ' ' X + i 
The expression (jlSp shows that the reflection coefficient tends to —1 as kod tends to 
0. This result is quite a striking one if one thinks of the extremely low concentration of 



material in this scattering experiment. The behavior of such a grid is equivalent to a perfect 
mirror at the vicinity of A = +00. For instance, for a/d = 1/1000 and for kod = 1/100, 
in normal incidence we find a theoretical reflection coefficient worthy of the best mirrors: 
R = |rp = 0.999739 ! Moreover the expression of the reflection coefficient r gives us the 
critical dimension of the radii of the wires. If a(A;o) is related to in such a way that 

M cos 6^0, d 

where F is some constant, i.e. 

d ttF 

a{ko) = — e kadcosOo 

then the grid, at the limit, does not behave neither as vacuum nor as a perfect mirror because 
the reflection coefficient is equal to r = Coming back now to the evanescent field and 

having in mind that kod ^ 1 we have /5„ ~ 2m\n\/d and therefore the expression written 
above ()13j) can still be simplifled: 



Eian - —7^ log (1 - e^^(^+^l^)l) E\Q) (18) 



III. DERIVATION OF THE CUT- WAVELENGTH 

In the preceding section, we have derived an explicit expression for the fleld diffracted 
by the wire grating, this allows us to derive the dressed Tq matrix j2(3| of a single grating 
layer. This matrix does not depend upon the characteristics of the incident wave, that is, 
the wavelength and the angle of incidence. 

Finally, in order to modelize the wire mesh photonic crystal, we derive the total dressed T 
matrix by adding a slab of air below the wire mesh grating. The matrix for a homogeneous 
slab of dielectric material with permittivity 1 between and —h is given by: 



cos(/5o/i) Pq^ sin(/5o/?-) 
-/?o sin(/5o/i) cos(/9o^) 
so that the total T-matrix for the grating and the slab is: 



(20) 



V 



^ ^ cos(/?o/i) + ^ sin(/5o/i) /3o ^ sin(/?o/i) . ^^i) 
—(3q sin(/5o/i) + J sin(/3o/i) cos(/5o/i) 



We have now characterized a basic layer of the photonic crystal. A general device is 
made of a stack of such layers. Due to the very weak evanescent fields, the transfer 
matrix of such a layered device is very well approximated by . Note that the influence 
of the evanescent waves can be easily derived from Eq. ^| 

In wire mesh photonic crystals, it has been well establish that there is a cut-wavelength, 
above which there is a band gap (this phenomenon is sometimes described as a plasmon 
frequency). Our point is to derive an implicit equation for the plasmon wavelength Ac. This 
wavelength corresponds to the edge of a band gap. A band gap is characterized by the fact 
that |tr(T)| > 2, when the wavelength is very large with respect to d, the transfer matrix 
is very near the identity matrix, consequently, the equation for the edge Ac of the last gap 
reads as: tr(T(A, 9)) = 2. This amounts to looking for /Jq = y- cos^^ which is solution of: 

cos{(3^oh) + -^M(3'oh) = l. (22) 

In order to solve this equation we denote:X = tan(/5Q/;,/2), and we obtain: 

1 2X 

whose solution is (3^L = X. Finally, we get that the plasmon frequency is given by the 
following implicit dispersion relation: 

/3^Ltan(/?^V2) = 1 • (24) 
Now, if we let = -rr^, this dimensionless number is solution of: 

' ATT ' 



27rx'=^tan(7rx^) = 1 , (25) 
that we have to solve numerically for any parameter ^. Three cases can therefore occur. 

• <^ 1 {i.e. h ^ for realistic size of rods, say a = 10~^d). In that case, we have 
~ 1/2 and, as a result, the cut-wavelength Aq is given by: 

X-^='^l^r^2hcos9. (26) 
Po 

which is nothing but the Bragg condition. As a conclusion Aq ^ (i (except for grazing 
incidence) which is, fortunatly, compatible with the homogenization process. 
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^ 1. In that case we can make equation 123 explicit by using the expansion 
tan(7rx'^) ~ ttx'^. We find: 



and we deduce an approximation Ag ^ of Aq 

Xf = nV2Lhcose (28) 

• ~ 1. We are therefore between the two previous cases. And in order to test the 
accuracy of our formula, we have to make some numerical experiments which is the 
topic of the next paragraph. 

In this paragraph, our aim is to evaluate the accuracy of Eq. |2S1 in some representative 
examples. Take the example of the stack of gratings each of them being made of thin circular 
metallic rods. In the following ^ = 0.005 has been taken for different ratios ^ and different 
kind of crystals defined by the angle ip in order to check the validity of our method. For 
this purpose, we plot the 0— order efficiency versus the normalized wavelength ^ in normal 
incidence as shown in Figs. 1-3. It is of importance to note that the cut- wavelength Aq 
derived from Eq. |^is extremely reliable even if | > 1 as shown in Fig. 1 for which f = 2. 
Conversely, the Bragg condition leads to absurd results except for very large ratios ^ = 2. 

Neglecting evanescent waves as we did to obtain the cut-wavelength, we can found an 
equivalent formulation as per: 



= 0, and 

y=0 



dE<^ 



= jE'^iO) (29) 



dy 

where [-jy^Q is the jump of the quantity within the brackets when crossing the plane y = 0. 
We have, indeed, E'^ solution of AE'^ + k'^E'^ = everywhere except for ?/ = and due to 
the invariance of the above relations with respect to x, we can look for solutions in the form 
E'^{x,y) = e^"o^t/(?/). The out going waves lead therefore to: 

Uiy) = \ ° (30) 
[ toe"'*^ ,y <0 

Eventually, by making use of Eqs. EHl we find expressions which are consistent with Eq. ITHt 
namely: 

^0 = 1 h-T to = l + ro . (31) 

i — tlJoL 

7 



As a conclusion, E'^ is solution in V 



(sense of distributions) of: 




y=o 1 



(32) 



and more generally for a stack of gratings 



7V-1 



AE'^ + klE' 



k=0 



y=—kh • 



(33) 



IV. CONCLUDING REMARKS 

In this paper, we derive a new expression for perfectly metallic gratings at long wave- 
lengths which leads to a very simple dressed matrix associated with one grating. In a second 
step we find a very accurate formula which gives the cut- wavelength of a very sparse photonic 
crystal even for cut-wavelengths of the same magnitude of the period of aforementionned 
crystals. This method can be easily generalized for non circular rods ; the characteris- 
tic length L cannot be derived in a closed formula and therefore necessitates a numerical 
computation. Moreover, the involved structures being infinitly conducting and completely 
transparent for the p— polarization, we think that it is possible to derive explicit formulae 
for three-dimensional structures. 
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Figures captions 



Figures captions 

Figure 1: Grid composed of thin infinitly metallic circular rods. 

Figure 2: 0— order efficiency Rq versus the wavelength A for a stack of 10 gratings in 
normal incidence for d = 1, a = 0.005 and h = 0.5. X^"^^^^ = 2h = 1 represents the Bragg 
wavelenghth, whereas Aq^ represents the ffist approximation given by Eq. |2Hland Ag is the 
cut- wavelength derived from Eq. 

Figure 3: 0— order efficiency Rq versus the wavelength A for a stack of 10 gratings in 
normal incidence for d = 1, a = 0.005 and h = 1. Ag"^^^^ = 2h = 2 represents the Bragg 
wavelenghth, whereas Aq^ represents the ffist approximation given by Eq. EE and Aq is the 
cut-wavelength derived from Eq. I2H 

Figure 4: 0— order efficiency i?o versus the wavelength A for a stack of 10 gratings in 
normal incidence for d = 1, a = 0.005 and h = 1. Ag'^'^^^ = 2h = 10 represents the Bragg 
wavelenghth, whereas Aq^ represents the ffist approximation given by Eq. |2Hland Ag is the 
cut-wavelength derived from Eq. |^ Note the widths of gaps for this very sparse structure 
(Filling ratio=1.5710~^!) 
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